Quantum inductance and negative electrochemical capacitance at finite frequency 
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We report on theoretical investigations of frequency dependent quantum capacitance. It is found 
that at finite frequency a quantum capacitor can be characterized by a classical RLC circuit with 
three parameters: a static electrochemical capacitance, a charge relaxation resistance, and a quan- 
tum inductance. The quantum inductance is proportional to the characteristic time scale of electron 
dynamics and due to its existence, the time dependent current can accumulate a phase delay and 
becomes lagging behind the applied ac voltage, leading to a negative effective capacitance. 
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Understanding dynamic conductance of quantum co- 
herent conductors is a very important problem in nano- 
electronics theory. When two quantum coherent conduc- 
tors form a double plate "quantum capacitor" , its dy- 
namic conductance G(lo) is given by the frequency de- 
pendent electrochemical capacitance^^ C^iuj), G(u) — 
—iujC^w), here u> is the frequency. At low frequency, 
C M (w) can be expanded in frequency and at the linear 
order, it is described 2 - by an equivalent classical circuit 
consisting of a static capacitor C M in series with a "charge 
relaxation resistor" R q . For a conductor having a single 
spin-resolved transmission channel, R q was predicted 2 
to be half the resistance quantum, R q = 1/2G where 
G = h/e 2 . The factor 1/2 in R q is of quantum origin 2 -^, 
and has recently been confirmed experimentally 5 . In the 
experiment of Gabelli et.al^, a submicron 2DEG quan- 
tum dot (QD) is capacitively coupled to a gold plate 
forming a double plate capacitor, where the QD con- 
nects to the outside reservoir by a single channel quantum 
point contact (QPC). The dynamic conductance G(cu) is 
then measured at 1.2GHz, and the data is well fit to 
the equivalent circuit characterized by two parameters 

(C^Rg). 

The experiment of Gabelli et.al^ opened the door for 
elucidating important and interesting physics of high 
frequency quantum transport in meso- and nano-scale 
devices. An important question is what happens to 
electrochemical capacitance at higher frequencies be- 
yond the linear uj regime, and in particular whether the 
two-parameter (C M ,i?q) equivalent circuit is adequate at 
higher frequencies to describe a quantum capacitor. It is 
the purpose of this paper to address these issues. 

In the following, we report a microscopic theory for 
high frequency quantum transport in a two-plate quan- 
tum capacitor. Our results show that to characterize its 
high frequency dynamic response, one needs — in addi- 
tion to C^, R q , a new quantity L q having the dimension 
of inductance. L q is found to have purely quantum ori- 
gin and will be named "quantum inductance" . Therefore, 
the frequency dependent electrochemical capacitance of a 
quantum capacitor C^(uj) is equivalent to a classical RLC 
circuit characterized by three parameters (C^, R q , L q ), at 



high frequency. Due to L q , electrons dwell in the neigh- 
borhood of the capacitor plates causing a phase delay. 
At low frequencies, the dynamic response is capacitive- 
like and voltage lags current. At larger frequencies when 
uj > 1/ yJC^Lq, inductive behavior dominates and volt- 
age leads current: in this case the quantum capacitor 
gives a negative capacitance value. It is, indeed, sur- 
prising that a quantum capacitor can give an inductive 
dynamic response. For the experimental setup of Ref[H, 
we estimate that when uj ~ 3GHz, the predicted high 
frequency effects should be observable. 

Let's first work out a simple expression for the fre- 
quency dependent electrochemical capacitance C M (w) fol- 
lowing the work of Buttiker— . We consider a two-plate 
capacitor similar to the experiment of Gabelli ef.a/.— : a 
QD — labelled I, and a large metallic electrode labelled 
II. Each plate is connected to the outside world through 
its lead and a time dependent bias vi,% is applied across 
the two leads. We consider small amplitudes of v% t 2 so as 
to focus on the linear bias regime. Under the action of 
such a bias, the two capacitor plates develop their own 
frequency dependent electric potential Utji(u>). The 
charge on plate-I is equal to the sum of the injected 
charge and induced charge: Qi = Q 1 ™ 3 + Q} nd . In the 
linear regime, the injected charge is proportional to bias 
Vx(u>): Q™ — e 2 Di{uj)vi(oj) where Dj(uj) is the general- 
ized global density of states (DOS) of plate-I at frequency 
uj. The induced charge, on the other hand, is in general 
related to a nonlocal Lindhard function 2 , whose calcula- 
tion is simplified by applying the Thomas-Fermi approx- 
imation. Within this approximation, the induced charge 
is proportional to the induced potential f/j on the plate: 
Qind _ —^jjj^^jjj where the minus sign indicates that 
the charge is induced. Putting things together, we ob- 
tain Qj = e 2 Di(uj)(vi — {//). Clearly, the same charge 
Qi can be calculated by the usual electrostatic geometric 
capacitance C : Qj = C (Ui — Uh). We therefore obtain 
a relationship: C (Ui — Un) — e 2 Dj(uj)(vi — U[). Ap- 
plying the same argument to plate-II, we similarly obtain 
—C (Ui — Un) = e 2 Du(uj)(v2 — Un). Finally, the same 
charge Qi can also be obtained from the definition of 
the electrochemical capacitance: Qi — C A1 (w)(ui — V2). 
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These three relations allow one to derive the following 
expression for C^(lo): 
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C„(w) C Dj{u) D u (cj) 



(1) 



This result resembles the one obtained by Biittiker for 
the static capacitance 2 . An important difference is that 
the frequency dependent electrochemical capacitance in 
Eq.([T]) is a complex quantity: its real part is a measure 
to the electrochemical capacitance and its imaginary part 
is proportional to the frequency dependent charge relax- 
ation resistance. 

To proceed further, we need to calculate the frequency 
dependent DOS D I:II {lu). Following Refl the gener- 
alized local DOS of plate-a of the capacitor can be ex- 
pressed in terms of Green's functions: 



dn a (u) 
dE 



dEf-f 
2tt huj 



[G r T a G a ] 
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where the subscript x labels space coordinates and 
D a (uj) = Tr[dn a (uj)/dE}. In Eq.©, / is the Fermi func- 
tion and / = f{E+) with E+ = E + huj; G r ' a = G r ^ a (E) 
is the retarded/ advanced Green's function at energy E 
and G r = G r (E + hu); T a is the linewidth function 
describing the coupling strength between plate-a with 
its lead. These quantities can be calculated in straight- 
forward manner when the Hamiltonian of the capacitor 
model is specified^. 

For the quantum capacitor of Refill, plate-I is a QD and 
plate-II is a metal gate. Since the metal gate has much 
greater DOS, i.e. Drj » Dj, we can safely neglect the 
Dji term in Eq.{T]). For a QD with one energy level Eq 
and connected to one lead, its Green's function G r — 
1/(E — Eq + iTl/2) where Tl is the linewidth function 
of the lead. With this Green's function, the frequency 
dependent DOS can be easily calculated from Eq.©, we 
obtain: 
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where A = AE 2 + Y 2 L /A, A± = (AE ± fuv) 2 + 
r|/4, and AE = E F - E . At resonance AE = 
0, we obtain Re(Dj(u>)) = [-x\n(4x 2 + 1) + 
2arctan(2a;)]/(27rr i x)/(a; 2 + 1) with x = huj/T L . Hence 
Re(Di(uj)) is positive for small x and negative for large 
x, i.e. there is a sign change. A similar behavior is also 
found for the system away from the resonance. Due to 
this sign change of Re(Dj), from Eq.© the frequency 
dependent electrochemical capacitance Cr = i?e[C M (u;)] 
can become negative. 

To be more specific, we fix the classical capacitance 
of QD Cq = IfF which is a typical value for QD with 
area of ~ 1/im 2 . Fig.l plots Cr = Re[C^(u))] (real part) 
and C'i = Im[C^(uj)\ (imaginary part) versus frequency 



for different values of Tl, by setting AE and tempera- 
ture to zero. We observe that Cr is positive at small 
frequency and becomes negative at larger frequency. For 
instance, Cr becomes negative at a "critical" frequency 
u! c ~ 10GHz for Tl = 10/ieV. This critical frequency 
can be smaller for smaller linewidth function T^. Wc 
note that it is not difficult to achieve Tl = 10/ieV exper- 
imentally: in Ref.8, Tl between lfieV to 5fieV has been 
realized. As will be discussed below, the effective in 
the experiment of RefH is tunable by a gate voltage so 
that the critical frequency at which the negative capaci- 
tance occurs can be even smaller. The inset of Fig.l also 
shows that as we increase u>, the imaginary part of C^^lo) 
starts from zero, reaches a peak value around oj c and 
then decays to zero. The negative capacitance at large 
frequency can be understood as follows. For a classical 
capacitor, a charge is accumulated across the capacitor 
induced by an external voltage. The current and volt- 
age has a fixed phase relationship: the voltage lags be- 
hind current with a phase tt/2. For a quantum capacitor 
at low frequency, there exists a charge relaxation resis- 
tance R q = h/(2e 2 ) for a single channel plate, therefore 
the charge build-up time is the RC-time trc = RqC^. 
For Cf, = IfF and R q = h/(2e 2 ), this RC-time is about 
trc — 13ps. If the external voltage reverses sign, the 
charge accumulation will follow the voltage and also re- 
verse sign in due time. When frequency is low, namely 
when uj < < 1 / trc — 77GHz , the charge build up follows 
the ac bias almost instantaneously just like a classical ca- 
pacitor, thus the capacitance is positive. When frequency 
is high, there is a phase difference between the ac bias 
and the charge build up. For frequencies comparable to 
1/trc, the charge build up can not follow the ac bias, 
thereby the capacitance may be negative. 

While the above argument explains why it is possible 
to have negative capacitance, it would indicate a critical 
frequency to be near 77GHz. Our results (Fig.l) show 
that the calculated uj c is actually much smaller. This is 
because there exists a second relevant time scale in the 
QD, i.e. the dwell time which is the time spent by elec- 
trons inside the QD. The dwell time Td can be calculated 
for specific systems 9 . Importantly, at resonance the elec- 
trons can dwell inside the quantum dot for a long time. 
For instance, for our QD with Tl = 10/ieV, we found 
r d = 260ps while t rc — 12ps (since C M = 0.9/F). In 
other words, » trc- Such a corresponds to a fre- 
quency of 4GHz, much less than 1/trc- In other words, 
when ac frequency is greater than 1/r^, the charges dwell 
inside the quantum dot and cannot respond to the ac 
voltage change. As a result, current becomes lagging be- 
hind voltage, leading to a negative capacitance. This pic- 
ture agrees very well with the numerical results (Fig.l). 

Having determined the general behavior of C^iw) for 
the quantum capacitor, in the following we determine 
how to simulate this quantity using a classical circuit. 
Expanding Eq. |T]) into a Taylor series to second order in 
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uj with the help of Eq.© at resonance, we obtain: 

cm - o,^-^^^ (4) 

where C M = C p (0) on the right hand side is the static 
electrochemical capacitance. This result is equivalent to 
that of a classical RLC circuit, as follows. For a classi- 
cal RLC circuit with capacitance C M , resistance R q and 
inductance L g , the dynamic conductance is 

G(uj) = -iwC„/(l - uj 2 L q C^ - iujCM . (5) 

Expanding this expression in power series of uj, we obtain 

G{uj) = -iujC^ + Lo 2 ClR q + iu?ClR\ - iw 3 C%L q . (6) 

Because for a capacitor G(uj) = —iujC^uj), we obtain the 
result that our quantum capacitance Eq. © is equivalent 
to the classical RLC circuit model of Eq.©. Comparing 
these two equations we readily identify R q — h/(2e 2 ) — a 
result first obtained by Biittikcr and co-workers 2 . Im- 
portantly, a new quantity — the equivalent inductance, is 
identified as L q = h 2 /(12tt c 2 Tl). In terms of dwell time 
Td and charge relaxation resistance R q , we obtain: 



L q = R q r d /12 



(7) 



where Td — AH/Tl. 

What is the reason that a quantum capacitor at fi- 
nite frequency needs to be modeled by a classical RLC 
circuit (instead of a RC circuit)? This is due to the 
role played by the large dwell time r d of the QD. When 
electrons dwells a long time Td inside the QD, the inter- 
action between electrons become an important piece of 
physics which, in our theory, is modeled by the induced 
self-consistent potential Uui discussed above. Such an 
interaction gives rise to the physics of induction, and re- 
sulting to the quantity L q of Eq.®. Indeed, the explicit 
dependence on Td by L q also confirms the important role 
played by the dwell time. Because L q is determined by Td 
as well as fundamental constants h and e, it is of purely 
quantum origin and can be called quantum inductance. 

Fig. 2 compares the fitting of classical RLC circuit 
Eq.©, with the full quantum result of Eq.©. They 
compare very well for the entire range of the frequency — 
if we treat R q as a function of uj. Indeed, while R q has 
so far been a constant h/(2e 2 ) as identified through the 
Taylor expanded Eqs. (|4l6[) , it is actually a function of uj 
by the more general expression Eq.©. The inset of Fig. 3 
plots the general R q — R q {uj) obtained numerically, and 
we observe it to be a slowly increasing function of a;. As 
expected, in the small frequency limit, R q {uj) recovers 
the result of half resistance quantum. For = 50(ieV, 
R q (uj) deviates from h/2e 2 at about 5GHz. We have also 
attempted using three constant parameters C M , L q and 
R q = h/(2e 2 ) into Eq.© to compare with the full quan- 
tum result of Eq. © , a reasonable agreement is obtained 
(inset of Fig. 2) although not as good as that shown in 
Fig.2. 



The situation is somewhat different for quantum in- 
ductance L q when the system is off resonance (AE 7^ 
in Eq.©). In this case the dwell time Td becomes too 
small to be relevant and another time scale becomes im- 
portant, namely the tunneling time r t for electrons to go 
in/out of the QD. The further away E is from Eq, the 
longer is t±. Hence in Eq.© Td should be replaced by Tt 
for off resonance. Our result shows that the fitting of full 
quantum capacitance C^(uj) using classical parameters 
Cfj,, L q and R q (uj) is still quite good for off resonance. 
This further supports the conclusion that the frequency 
dependent quantum capacitance can be described by a 
classical RLC circuit with static electrochemical capac- 
itance, charge relaxation resistance, and a quantum in- 
ductance. 

Finally, we perform a numerical calculation of the dy- 
namic conductance for the device structure of Ref|5|. In 
terms of scattering matrix, the DOS of Eq.© for a ca- 
pacitor can be re-written as^£ 



dEf-f 
2tt H 2 uj 2 



(E+)sll(E)] (8) 



witb£ s[ L {E) = (r - e^)/(l - re 4 *), tj> = 2irE/A, 
r 2 = 1-Tqpc and T QPC = 1/[1 + exp(-(Vg + V )/AV )\ 
which is the transmission coefficient of the QPC in the 
experimental setup^. Fig. 3 shows the dynamic conduc- 
tance G{oj) vs gate voltage V g using our theory presented 
above. When uj = 1.2GHz (open circle), our results 
agrees very welt^ with the experimental data of Ref[H. 
When uj = 3GHz (open square), our theory predicts that 
the imaginary part of G(u>) which is the electrochemical 
capacitance, goes to negative. For even larger frequency 
uj = 5GHz the effect is more significant. To understand 
why one can observe a negative capacitance at small fre- 
quency such as 3GHz, we note that since electron enter- 
ing the QD has to first pass the QPC—: this QPC serves 
as a barrier (with an effective barrier height 1/T) that 
is controlled by the gate voltage. At small gate voltage, 
Tqpc is nearly zero and goes to one at large V g . Hence 
the effective T for small gate voltage is quite large making 
uj c much smaller. Since the experiment of Ref[^ is per- 
formed at uj — 1.2GHz, we assume that it is not too dif- 
ficult to push the frequency to 3GHz so that the effect of 
quantum inductance can be observed experimentally In- 
deed, a single-wall carbon nanotube transistor operated 
at 2.6GHz has been demonstrated^ and measurement of 
current fluctuation at frequency from 5 to 90 GHz has 
been reported^. 

In summary, we found that at finite frequency, a quan- 
tum capacitor consisting of a quantum dot and a large 
metal conductor is equivalent to a classical RLC circuit 
with three basic parameters: a static electrochemical ca- 
pacitance C^, a charge relaxation resistance R q , and a 
quantum inductance L q . It is found that L q ~ R q T where 
r is the characteristic time scale for the quantum dot such 
as the dwell time T d or the tunneling time r t . Due to the 
phase delay by the quantum inductance, the dynamic 
current can lag behind of the applied ac voltage, giv- 
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ing rise to a negative capacitance. Our numerical results 
show that this effect should be detectable experimentally 
using the present device technology. 
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FIG. 1: Frequency dependent capacitance C M (w) versus fre- 
quency. Main figure is for Cr and inset is for Cj. Here 
Tl — WfieV (for the curve with the open circle), 50/ieF 
(open square), 100/j.eV (open triangle), and 200/ieV (open 
diamond). Here C = IfF. 



FIG. 2: Comparison of the full quantum capacitance C M (w) 
(solid line) to that obtained by the classical RLC circuit model 
(dotted line) for Tl = 50fieV. Here Re[C^{u:)] is indicated 
by open circle. Inset: similar fit but using three constant 
parameters. 
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FIG. 3: The dynamic conductance G(u>) (unit e 2 /h) versus 
gate voltage at different frequencies uj = 1.2,3,5GHz. Solid 
line — Re[G]/uj, dotted line — Im[G]/ui. Other parameters: 
A = 500mA", a = 3000, V = -0.85, AVb = 0.003, C = 
4/F, and temperature 50mK. For purpose of illustration, 
we divided G(u) by u). Inset: frequency dependent R q (unit 
h/e 2 ) vs frequency. 
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